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This report presents the results of a study whose principal 
objective was the development of a mathematical model of a 
hingeless rotor for a tilt-rotor aircraft by synthesizing a 
set of wind tunnel data obtained on a model of the rotor. A 
secondary objective was to incorporate this rotor math model 
into a real-time flight simulation model of a hingeless rotor 
XV-15 tilt rotor aircraft emd determine the effects on the 
aircraft trim and stability as compared to a former model de- 
veloped under a previous contract. 

The wind tunnel data used in the synthesis was obtained on a 
1/4.622 scale model of the Boeing Model 222 tilt rotor air- 
craft tested under an earlier phase of the contract. The test 
generated sufficient data to define the rotor behavior over 
the range of flight speeds anticipated for the modified XV-15 
aircraft. 

The study was aimed at developing a set of equations that 
(1) would represent the rotor behavior as determined by the 
wind tunnel test data^ (2) could be evaluated rapidly by a 
computer so that the computation cycle time requirements of 
real-time piloted simulation would be met. 

Two approaches to developing the math model equations were 
tried before the final method was selected. The first 
approach was to apply the techniques of statistical linear 
regression to the wind tunnel model rotor data. The second 
approach was to develop a simplified analytical model for the 
hingeless rotor/ correlate the predictions of the analytical 
model with the measured data and then use the correlation 
functions with the analysis to yield corrected values. These 
approaches are presented as appendices. 

The actual method adopted was to apply a systematic curve- 
fitting procedure. The resulting equations together with an 
extensive set of graphs of the rotor derivatives as obtained 
from the test data are presented. Numerous plots of the math 
model results against the corresponding wind tunnel data are 
provided. These show that the math model equations reproduce 
the test data very well. 


The equations for the rotor were programmed into a previously 
developed flight simulation math model of the hinge less rotor 
XV-15 and the effect on aircraft trim and stability observed. 

It was found that trim is about the same as that calculated 
with the former rotor math model. The aircraft control deriva- 
tives are sufficiently different, however, that control phasing 
and rescheduling are required in order to maximize control 
efficiency. 
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1.0 INTRODUCTION 

In many applications of rotor technology^ there arises a need 
to be able to compute rotor forces and moments rapidly and 
accurately. One such application is real-time flight simula- 
tion where the rotor effects must be updated at intervals of 
less than €0 milliseconds. Existing finite element programs 
or even modal programs which solve the blade dynamic equations 
cannot satisfy the requirement for fast computation. In 
addition, the accuracy of the answers is limited by the so- 
phistication of the analytical model, greater accuracy de- 
mands more detailed modeling which increases computation time. 

Data Bank 


One solution that has been used extensively in the past is to 
construct a bank of experimental or, less desirably, analytically 
derived data. Values of rotor hub forces and moments are 
then obtained by a rapid search and Interpolation procedure at 
specified values of the rotor flight parameters. While this 
approach may be faster than a purely analytical computation, 
the range of flight parameters over which the Interpolation 
must take place (airspeed, collective, shaft angle of attack, 
tip speed, cyclic control) and the need, in some cases, for 
quadratic rather than linear interpolation, place large de- 
mands on data storage and retrieval. 

Functional Representation 

A more promising approach is to synthesize the data obtained 
from test or analysis into a set of equations for each vari- 
able. For example, it may be possible to develop an equation 
for rotor normal force Involving basic parameters such as y, 

®.75» Of Ai, Bi. The resulting equations may be lengthy and 
complex, involving products and trigonometric functions of 
the basic variables, nevertheless, they can provide a very 
rapid means of calculating the forces or moments. The data 
synthesis approach may also yield benefits beyond the rapid 
reproduction of the data, in that functional relationships of 
a general character may be recognized and may be useful in 
understanding the physical behavior of the rotor. A further 
advantage of having a set of equations to represent the rotcr 
forces and moments is that the equations may be used to esti- 
mate the behavior of a similar but untested rotor, and provide 
guidance in the preparation of a cost-effective wind tunnel 
test plan. In fact, once the general form of the equations is 
established, the wind tunnel testing could be regarded as a 
calibration activity. 
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The functional synthesis approach was adopted during a recent 
study (Reference 1) which involved development of a mathemat- 
ical model of the XV-15 aircraft with a hingeless rotor for 
piloted simulation. In this study/ wind tunnel data on a 
full-scale hingeless rotor obtained under a previous phase 
of this contract (Reference 2) was reduced to a relatively 
simple set of equations. Gaps in the wind tunnel data re- 
sulting from pressure of time and limitations on the tunnel 
speed were filled by calculation. 

Satisfactory experience with this functional model of the 
rotor led to the conclusion that it would be worthwhile pro- 
ceeding to a more definitive and complete model based on the 
extensive data/ reported in Reference 3/ that was acquired by 
testing a 1/4.622 scale model of the same rotor. The develop- 
ment of this new functional representation and the definition 
of a broad general procedure for performing functional syn- 
thesis of rotor test data is the subject of the present report. 

Present Studies 

Three approaches were explored. The first was an attempt to 
apply multivariable linear regression techniques to the data. 
This was motivated by the need to process the large amount of 
data available in the shortest time. The attempt met with 
limited success and is described in Appendix A. One of the 
difficulties encountered using the regression method was in 
making the correct choice of regression variables. From the 
previous work it was known that functions like ycoso/ 
p^Cfsina/ etc. would be required. There are/ of course/ many 
possible combinations of functions of the basic variables 
that will yield a satisfactory correlation/ and a large fac- 
tor in the success of the regression technique is making the 
correct choice of functions. It waS/ therefore/ concluded 
that theoretical guidance was required. A simplified theo- 
retical analysis was then written to identify these functions. 
Later/ the analysis was tried in a predictor-corrector tech- 
nique wherein the theoretical predictions were to be corre- 
lated with the test data and a correlation function generated. 
This approach was not pursued when it became evident that the 
analysis did not correctly predict trends with advance ratio. 
Results of the analysis and the approach are included as 
Appendix B. 

The approach finally adopted was to develop curve fit equa- 
tions for the derivatives 3CHF/3a/ 3CNF/dAi/ etc. These were 
then used to reduce the data to a reference set of values of 
Cp/ Cnf> Csf# etc./ that were dependent on thrust coefficient/ 
advance ratio and angle of attack only. This dependence was 
then curve fitted. A full description of the approach to to- 
gether with extensive correlations with test data is presented 
in Section 3. 
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2.0 DATA BASE 

The data which forms the basis for the study reported herein 
was obtained during extensive tests of a X/4.622 Froude-scale 
model of the Boeing Model 222 tilt rotor aircraft. The test- 
ing and results are reported fully in the four volumes com- 
posing Reference 3. A brief description of the test and the 
nature and extent of the data is provided in this section. 

Figure 2.1 shows the tilt rotor model installed in the Boeing 
Vertol wind tunnel. The test was designed to provide force 
and moment data on both the airframe and rotor over the ex- 
pected normal range of flight speeds and attitudes. Data was 
obtained at seventeen points in the flight spectrxim as shown 
by Figure 2.2. At each of these points, variations were made 
in fuselage angle of attack, yaw angle, collective pitch, 
cyclic pitch, wing flap setting and rotor rpm. 

The extend of the data is indicated by Table I, which pre- 
sents a summary of the configurations tested. The selection 
of test points was made in such a way that a comprehensive 
set of data was obtained for all potential flight conditions 
from hover through transition to cruise flight at simulated 
speeds up to 300 knots. Rotor force and moment data on each 
rotor was obtained from nacelle-mounted balances. The blades 
were strain gauged to provide flap and chordwise bending mo- 
ment information. 

Data Inspection 

Before beginning a synthesis of the data, the plotted results 
contained in the four volumes of Reference 3 were inspected 
in order to identify possible bad data points attributable to 
such effects as zero shifts. This process resulted in a 
decision to utilize the data measured on the left-hand rotor 
only since the right-hand rotor data was questionable in some 
areas due to strain gauge failures. The acceptable portions 
of the data on the right rotor were used, however, to cross 
check the behavior of the left-hand rotor. 

Data points that were called into question were inspected 
further using the microfiched records of the balance force 
data. As a result of this, the data points were either re- 
tained or discarded. The good data points were then punched 
out on cards and read into a computer program for formating 
and analysis during the next processing step. 
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3.0 REDUCTION MSTHOD AND COMPARISON WITH TEST DATA 

Rotor forces and moments (referred to axes fixed in the rotor] 
depend on seven independent parameters - collective pitch 
(9.75) » lateral cyclic (Ai) , longitudinal cyclic (Bi) , shaft 
angle of attack (a) « rotor sideslip angle rotor rpm and 

airspeed. For data reduction purposes it was assumed that 
any of the force or moment coefficients cculd be represented 
by an equation of the form 


(Ai-A^ 


(Bi-Bi 


+ 3Cp/3t(» ^ * 3Cp/3RPM (RPM-RPM^^^) 


( 1 ) 


in this equation, CpR is the value of the coefficient at 
reference values of cyclic, angle of attack and rpm. The 
reference values are the nominal values at which the rotor 
was first trimmed to minimum blade loads and then collective 
pitch varied. Thus the coefficient Cp|^ may be written 

^'®®*“ref' ^1 ref' ®1 ref' *^”ref^ 

The variation of these nominal values with the effective rotor 
advance ratio, ucosa, are presented in Figures 3.1 and 3.2i 
Deviations from these reference conditions are then accounted 
for by the derivative terms. 

In order to establish the derivatives SCp/SAj^, etc, a run-by- j 

run, data point~by data point inspection of the plots of I 

Reference 3 was conducted. It' was found that quite often j- 

shifts occurred in one or more of the parameters that were to I 

be held constant. For example. Run 41 was a run in which | 

lateral cyclic was varied with the rotor at 80* angle of 
attack and 80 knots airspeed. For the first six data points, 
the longitudinal cyclic (B^) was held fixed at 5.0 degrees. 

However for data points 7 through 11, the value of Bj^ jumped 
by 0.3 degrees. By carefully noting such shifts and by using 
only those data points that met the criteria for a derivative, 
it was found that a reasonably consistent set of derivatives 
could be obtained. The values of the derivatives were plotted 
against ucosa and curves hand-faired through the points. 

These curves were then fitted by a polynomial in ucosa of the 
form 



3CF 

3X 



(ucosa) 


i 


(3) 
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Iftl 


where X represents a, and rpn. The values of these 

coefficients are listed in Tables II through IX. The varia- 
tion of the derivatives with vcosa are presented in Figure's 
3.3 through 3.38. 

3.1 Determination of the Values of the Coefficients/ 


The equations for the derivatives/ together with the reference 
values of cyclic/ angle of attack and rpm obtained from 
Figures 3.1 and 3.2 were then used in Eguati.'.n (1) to compute 
values for CpR for those runs in which collective alone was 
varied. The expectation was that the Cpi^ could be a linear 
function of e.75 or Cxr« Figures 3.39 through 3.44 present 
the results of this procedure and it can be seen that the 
expectation of linearity is confirmed. 

The next step in the reduction process was to obtain a set of 
equations that would reproduce the behavior of Cpp vs at 
the different flight conditions/ i.e. since 

^^FR 

C * C 4- ^ C ' (4) 

'"FR '“FO ^TR 

the dependence of Cpo dCpR/SCxr on w must be established. 

The values of Cpo and dCpp/dCxp were plotted against wsina 
since previous work indicated that the dependence would be of 
the form 

CpQ « f (usina,u) 


" ■ f (ysina/w) 

^^TR 

Figure 3.45 shows the result of this for the case of pitching 
moment at zero thrust/ Cpi^Q. It can be seen that there is 
insufficient data at fixed u and different values of a to be 
able to establish the form of the functions. Various combina- 
tions of vsinof ucosa were tried in an attempt to obtain 
smooth variations of the Cpg and 3Cpp/aCxp that would include 
all the data points. None was successful. The most satis- 
factory approach was to plot the data against u/ fair a line 
through those points that lay closest to the nominal shaft 
angle of attack schedule (Figure 3.1)/ and then fit an equa- 
tion to the line. A typical result is shown in Figure 3.46. 
This process was repeated for all the rotor force and moment 
components . 


I 

[ 

[ 


r 


} 

i ' 
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3 . 2 Ma th Model Equations 

The final equations for the rotor forces and moments are 
lengthy since they are combinations of polynomials in v and 
ucoso whose coefficients change depending on which portion of 
the range of p or pcosa the rotor is operating in. For this 
reason, the entire set of equations will not be presented 
here. The equations are given in the form of a small computer 
subroutine in Appendix C. 

The general equation used to calculate the forces and moments 
is Equation (1) viz 


C 


P 


where 




- =PR ^ I57 '^-*1 ref> * <*r®l ref> 


+ 3Cp/3i|» ^ + 3Cp/3RPM (RPM-RPM^^^) 

3C„„ 


'NF 


'SF. 


'PM, 


TR “ 



'PR “ 


'Sr 


+ 

'Sfr 


... 0 °''TR 

^^SFR 

C + ^^pmr 

^^0 ^Sr 


'TR 


'TR 


'YM, 


3C 


= C 


YMR 


'Sr 


'TR 


3CF 

3a 


(a-o^ref 


) 


(5) 

( 6 ) 

(7) 

( 8 ) 
(9) 

( 10 ) 


in which the quantities 3 C.tij^/3 0 ^ 75 , 0 q» CpQ, 3Cpp/3Cmp, etc, 
are obtained from piecewise curve fit equations with p or 
pcosa as argument. The derivatives, 3C /3A. , ac /3B , 
3C_„/3a , etc., are obtained from 1 nf i 

rM 

where the aj^ are listed in Tables II through IX, and are 
furnished by the array X(I,J,K) in the subroutine. 
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The equations that define the reference quantities oref» 

^Iref' ®lref ^^^ref' Presented in Figures 3-1 and 3-2, ares 

“ref “ " 322.82 + 770.412^ - 1446.342^ + 1449.92^ 


- 545.1152- 


where 


2 ° ucosa 

*lref ■ ®lref ‘ 5-'>' 


and the reference rpm is a function of nacelle angle, 


RPMref = 1185 


» 1185 - 7.89(45-ijj) 




The calculation procedure is illustrated by the following 
example . 

If '^<ycoso< . 015 , 


» (11.9332 + 0.882)/1000 


0Q = -1.5 + 106.6672 

and equation (5) yields the reference thrust coefficient C . 
Equation (11) is now used to calculate the derivatives and^^ 
equations (12) , (13) and (14) used to compute the reference 
quantities. The total thrust coefficient is then calculated 
for the given values of cyclic, angle of attack, yaw angle, and 
rpm. 

For the same range of ycosa the power coefficient at zero 
thrust, CpQ is 

Con = (-035 - 9.667Z)/1000 


.094 - .2672 
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The total power coefficient is then obtained in the same way 
as thrust. 


If 0<u<.51 the equation for normal force coefficient at zero 
thrust, C^^pQ, is 

1000 X Cj^p = -1.45 + 1.1164U + 46.159u^ + 46.153u^ 


-4118. 75u^ + 27714. -74786u® + 89304u^ 


8 


-39318U 

and the gradient with thrust coefficient is given by 


dC 

dC 


NFR , 
TR 


-.0087 + 2.122u -10.95w^ -34.52y^ + 420u^ 
-1043u^ + 520. 8u^ + 1029u^ -966u® 


Using these quantities in equation (7) the reference normal 
force coefficient is calculated and finally the total normal 
force obtained from equation (1) . The Scune procedures are 
used to calculate sideforce, pitching moment and yawing 
moment . 


3.3 Comparison with the Test Data 

The equations described in 3.2 were used to generate estimated 
values for the left rotor forces and moments at all the condi- 
tions investigated during the entire wind tunnel test. An 
extensive series of plots of the estimated and actual forces 
and moments were then made . 

3.3.1 Collective Sweeps 

The first test of the ability of the equations to reproduce 
the test data is to see how well the forces obtained during 
collective sweeps are estimated, since this is the data which 
was the most difficult to collapse. Figures 3.47 through 
3.52 present comparisons of estimated and test values of 
Cf^, Cp, Cfjjf-, Cgp, Cpjj and Cy^f respectively. Thrust is esti- 
mated to within 10% of the measured thrust at the same value 
of collective pitch. The estimated variation of power coeffi- 
cient with thrust coefficient is also satisfactory at all the 
test conditions with the exception of hover at high values of 
Ct(. 013) where the estimate is about 9% too low. Normal force 
estimates are also acceptable bearing in mind that, unlike 
power and thrust, normal force is very sensitive to cyclic 
and the accuracy achieved depends on how well the values of 
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the cyclic derivatives are fitted. Similar remarks apply to 
the comparisons of estimated and actual values for sideforce, 
pitching moment and yawing moment. 

3.3.2 Detailed Comparisons at Points in the Transition 
Corridor 

Figures 3.53 through 3.135 present a detailed comparison of 
the estimates obtained from the math model equations with the 
test values at three points in a typical tilt rotor transition 
corridor: 

1 . Hover 

2. 45 knots and 90 degrees nacelle angle 

3. 180 knots and 15 degrees nacelle angle 

The comparisons are presented in the form of copies of the 
original wind tunnel data plots with the estimated values 
superimposed. 

In hover. Figures 3.53 to 3.70, the rotor response to lateral 
and longitudinal cyclic is estimated very well considering 
the degree of scatter that exists in the measured data. The 
behavior with collective is reproduced nearly exactly. 

For the early transition test condition of. 45 knots and 90 
degrees nacelle angle. Figures 3.71 to 3.106, the estimated 
variation of the forces end moments with lateral cyclic, A]_, 
is generally in good agreement while the response to longi- 
tudinal cyclic, Bi, is in very good agreement with the ex- 
ception that the estimated yawing moment is displaced from 
the measured values. The correct slope of yawing moment with 
Bi is predicted correctly, however. The dependence of the 
forces and moments on 675 is estimated well, again with the 
observation that yawing moment is displaced. For those runs 
where shaft angle was varied, Figures 3.89 through 3.94, it 
can be seen that the math model shows a weak dependence on 
shaft angle, o. This is attributed to the use of the parameter 
vicosa in the equations since ycoso changes rapidly near a = 

90 degrees. Comparison of the estimated and test behavior 
of the data during yaw angle sweeps shows that the effect of 
yaw angle is not correctly estimated for normal force and 
pitching moment. This defect is probably due to the derivative 
curve fit in this region. The effect of rpm changes is 
reproduced well. 

Figures 3.107 through 3.135 present the comparisons for a l\igh 
speed end of transition condition, 180 knots and 15 degrees 
nacelle angle. Response to longitudinal cyclic is estimated 
extremely well for thrust, normal force, sideforce, pitching 
moment and yawing moment. The response of power coefficient 
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to Bi is too exaggerated, however. The angle of attack 
sensitivities are reproduced very well and predicted yaw 
angle effects are also satisfactory. The effect of rpm vari- 
ations are less satisfactory for this condition. 

From the comparisons presented here it is concluded that the 
curve fit equations reproduce the test data sufficiently 
accurately for simulation purposes. The incorporation of the 
equations into the simulation math model of Reference 1 and 
the results are discussed in the following section. 
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4.0 EFFECT ON AIRCRAFT TRIM AND STABILITY 

The math model equations for the rotor described in Section 

3.0 were incorporated in the simulation model of the hinge- 
less rotor XV-15 aircraft. This simulation model is detailed 
in Reference 1. The new equations, in the form of the sub- 
routine given in Appendix C, replace the set of equations 
described in that Reference. An exception is that the pitch 
and yaw rate derivatives were retained since the test data on 
which the new equations are based does not contain rate effects. 

The new rotor equations yield the forces and moments in a set 
of axes fixed in the rotor rather than the wind axes system 
used in the former equations. Some existing transformations 
from wind to body axes were therefore eliminated and others 
added. 

4 . 1 Comparison with the Former Rotor Model 

The simulation math model incorporating the new rotor equations 
was used to assess the effect of the updated rotor representa- 
tion on the trim and stability of the XV-15/hingeless rotor 
aircraft and to ascertain whether control schedule changes 
would be required. Tables x through XIV present comparisons 
of aircraft trim attitude and control settings required and 
derivatives at selected points in the transition corridor. 

From Table X it can be seen that the introduction of the up- 
dated rotor representation does not substantially change the 
values of trim pitch attitude ( 3 ) , stick position (6g) or 
throttle setting values of lateral (Aj^) and 

longitudinal (Bj^) cyclic (referred to the classical wind 
axes) are changed, however.* With the new rotor math model, 
approximately twice as much lateral cyclic is required and 
about ha lf as much longitudinal cyclic. The total cyclic, 

^ is about the same with both rotor representations. 

The different cyclic requirements reflect the different sensi- 
tivities to cyclic and angle of attack given by the two models. 
Since the new rotor representation has been shown to be in 
very good agreement with the test data, a re-phasing of the 
cyclic inputs to maximize control force and moment is 
indicated. 

Tables XI through XIV present preliminary comparisons of the 
aircraft derivatives obtained using the present and former 
rotor models. Because of the altered cyclic sensitivities 
some of the control derivatives are changed and the control 
system gains and schedules require reworking to provide 
acceptable flying qualities. The controls-fixed derivatives 
M_, N„, X„, y„, Z.,, for the 60 knot, i„ » 75* condition 

p 4 r U V W N 
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are comparable with the two models. The off-diagonal deriva- 
tive, Yd, side force due to roll rate, is changed in sign 
compared to the present model. The reason for this has not 
yet been established. For the 140 Kt, i)f » 60** case (Tables 
XI and XII) , the principal derivatives Lp and differ from 
the former model. The source of these differences is yet 
to be uncovered but a preliminary examination indicates that 
the new model exhibits a higher sensitivity of thrust to 
combinations of angle of attack and effective rpm changes 
such as occur during rolling, than did the former model. 
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5.0. CONCLUSIONS AND RECOMMENDATIONS 

The objective of the study was to develop a functional repre- 
sentation of a hingeless rotor based on extensive wind tunnel 
data and to incorporate this representation into a real-time 
simulation math model of a hingeless rotor XV-15. Changes in 
flying qualities caused by the new rotor representation were 
to be determined and noted. 

Conclusions 


1. A mathematical representation of the rotor data was devel- 
oped through a systematic curve-fitting approach. The 
results show that the curve-fit equations reproduce the 
wind tunnel data very well. 

2. The rotor equations were incorporated into the simulation 
model. Preliminary evaluation of trim and derivative 
data show that cyclic rephasing is required in order to 
optimize the aircraft response to control. Stick- fixed 
behavior of the aircraft is generally comparable to that 
obtained using the former rotor representation. 

3. Extensive plots of the rotor derivatives as functions 
of the effective advance ratio have been presented. 

These are of considerable value in themselves since they 
may be used to assess rotor behavior for similar configu- 
rations . 

4. A curve-fit representation of the rotor, although lengthy, 

satisfies the computational speed requirements of real- 
time simulation just as well as the equations in a 
functional representation. This is because the functional 
representation, although simpler in form, involves comput- 
er evaluations of transcendental functions, e.g., Tdn~‘x, 
sinrt et require as much time as the evaluation of 

polynomials in the curve-fit approach. 

5. An attempt was made to automate the process of data syn- 
thesis using multivariable linear regression. In princi- 
ple, if a sufficiently large number of candidate corre- 
lates are available to the regression algorithm then a 
successful synthesis should result. However, the likely 
correlates must be supplied to the computer and the selec- 
tion of the correct correlates requires a considerable 
amount of skill, if not luck. The final equations produc- 
ed by the regression, while they may fit the particular 
set of data used, are not necessarily of a general enough 
nature that they can be extended to other rotor configu- 
rations by simply altering constants in the equations. 
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Recommendations 


1. Probably the best overall approach to the problem of pro- 
viding a fast/ but general/ algorithm for rotor represen- 
tation in real-time simulation is to develop a simplified 
analysis whose predictions can be calibrated against wind 
tunnel data. These analytical predictions would then be 
operated on by the calibration functions to produce esti- 
mates of the rotor forces and moments. A step in this 
direction was made in this study but more work is 
required. 

2. A further recommendation is that during the preparation 
of wind tunnel test run schedules/ every effort should 

be made to select the runs so that the natural parameters 
of the rotor are systematically varied over a broad enough 
range that the data synthesis process is facilitated. 

For exeunple/ when testing at different rotor shaft angles 
at fixed airspeed, such as is required for transition 
corridor definition, it is desirable to do this at a ser- 
ies of fixed values of rotor rpm even though the aircraft 
will not necessarily operate at some of the rpm values. 
This will, however, ensure that data exists over a wide 
enough range of values of u, ysina, wcosa that curves may 
be confidently drawn through points of constant psinoi, 
for example, and equations fitted to them, 

3. The final recommendation arises from the difficulties 
experienced in extracting rotor derivatives from the 
test data. The wind tunnel run schedule was purposely 
designed to acquire rotor derivative data directly, by 
making runs in which one control at a time was varied 
while the others were fixed. 

Unfortunately, during many of the runs, shifts in cyclic, 
collective pitch or shaft angle occurred. This made the 
job of extracting rotor derivatives very laborious since 
each point had to be inspected to make sure that the 
fixed controls had actually remained fixed during the 
run. It is therefore recommended that in future tests of 
this nature every effort should be made to ensure that 
the control settings remain at the desired levels. This 
may require the use of more complex and more accurate 
rotor control systems than the simple position- feedback 
system used in this test. 
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TABLE I. CONFIGURATIONS TESTED (Concluded) 



Table II. Coefficients for Thrust Coefficient Derivatives 
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Table HI. Coefficients for Power Coefficient Derivatives 
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Table V . Coefficients for Sideforce Derivatives. 
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Table VIII . Coefficients for Chord Bending Moment Derivative 
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TABLE X. COMPARISON OF AIRCRAFT TRIM CONDITIONS AS COMPUTED 
WITH THE PRESENT AND FORMER ROTCR MODELS 




D210-11505-1 


c 



a^ 


n 

>0 


c 


00 

/ 


K u: 

r 





in 

a> 

in 


9m 

/ 

»“ 

Ll o 

V 


n 


o 

in 

ro 

9^ 



r nO 

c 


NJ 

cn 

cn 

o 

o 

c 

n 

C 

o 



c 



• 

• 

• 

• 

# 

• 

• 

• 

•/ 

• / 

• 

o o 




1 



cu 

• 

1 

/ 

• 





o 

9m 

o 

o 


O 

n 

o 

/oo 

/o 

o 


X 

o 

tn 


o 


O 

in 

o 

/o> 

/ o 

o 


V 

o 

•-4 

(n 

o 


c 



r ^ 

o 

c 

H M 

> 

o 

n 


c 

fv 

o 

o 


CVJ 

c 

w 

• O. 



• 

1 

« 


1 


/ 

1 ' 

• 



rv r> 

^ o 

0% ^ 

m o 


^ >o 

^ (U 

ro ^ 

o 


2/r / 

a i ^ 


n ♦ >c 

^ I <u <u o 

^ / o rv o 

^ / o o o 


rs4 

M 

o 


>6 

o 



/cu 

/ 

r - 

n 


o 

(n 

r> 

o 

00 


/in 

/ ^ 

OG 

V 

o 



o 


o 

/ a* 


o 

z 

o 

A 

fU 

m 

c 

m 

in 

• 



^ o 

• 

o 

• 


•/ ■ 


> 

> 

(U 

0^ 

00 

«<# 

n 

/n 

m 

in 

n 



•-4 

fU 

M 

00 

o 

9>4 

/<Vi 

/-e 

n 


in 

o 

V 


O 

in 

o 

o i 

r i 

r ^ 


o 

c 

c 

i: 


O 

o 

o 

O / 

^ / 

^ o 

o 

o 

o 

Q 


»• 

• 

' • 

• 

• / 


• 

• 

• 

• 

• 





• 

/ 

it/ 

• 



• 

• 




in 


X 






/ 



in 



* 


n 


X 

o 

o 

00 

o 

01 


o 

o 

o 

o 


0) 

w4 


09 

o 


fN 

o 

o 

/ 

>0 

o 


w 

w 





N 

c. 

(K 


c 

C J 

/ C 

(V 

c 


c 

c 

T ‘ 

t 

II 

n 


c 


c 

c 

0 y 

' c 


c 

««■ 

c 


* 

> 




• 

• 

• 

• 

• / 

• 

• 

• 

0 

• 

• 

* 

X 

# 

m 




^ • 

• 

w 

• 

1 

• 

• 

• 

• 



s 

iS 




% 


»/ 









0 

• 






/ 







* - 


o 

u 






/ 







** 








X 













a: 

Hi 

oe 

c 







c 









111 








• ” 






a 

u 

X 







z 

, 





o 

O 

o 


(L 

o 

ce 


> 


0m 


31 


TABLE XI. DERIVATIVES AT 140 KT, i„ = 60®, «„ » 0® , - FORMER ROTOR MODEL. 
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TABLE XII. DERIVATIVES AT 140 KT, i •= 60*, « = 0* - PRESENT ROTOR MODEL 
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TAILE XIII. DERIVATIVES AT 60 KT, I - 75*, - 40" - FORMER ROTOR MODEL 
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TABLE XIV. DERIVATIVES AT 60 KT, = 75®, «„ = 40*- PRESENT ROTOR MODEL 
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Figure 2.1. 1/4.622 Scale Model Installed in the Wind Tunnel 
Test Section 
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Figure 3.1. Variation of Reference Angle of Attack With 
li Cos 0( 
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Figure 3.2. Definition of Reference Values of Cyclic and RPM 





Figure 3-3. Variation of 3C with ptcosa. 
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Figure 3-4 • Variation of 3C— with uCOSa 
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Figure 3-8. Variation of 3Cp with ycosa 
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Figure 3-13. Variation of 3C„_ with poosa 
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Figure 3-16. Variation of with pcoso 








Figure 3-17. Variation of with vicosa 
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Figure 3-18. Variation of 3C„ with ycosa 
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Figure 3-22. Variation of 3C with iicosa 







Figure 3-23. Variation of with ycosa 
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Figure 3.25. Variation o£3Cb„/DRPM with ucosa 












Figure 3-28. Variation of with pcosa 
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Figure 3.30. Varation of 3 Cvm/ 3RPM with ycosa. 




Figure 3-31. Variation of 3FBM with pcosa 






Figure 3-32. Variation of 3FBM with pcosa. 






Figure 3-33. Variation of 3FBM with pcosa. 










Figure 3-25. Variation of 3CBM with ycosa. 
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Figure 3-37. Variation of 3CBM with ycosa 
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Figure 3.38. Variation of SCBM/3i|) with ucosa. 


U 






M80L 

RUf 

* 

27 

o 

43 

0 

51 



D2101 1505-1 



Figure 3.41. Corrected Normal Force Coefficient Versus 
Corrected Thrust Coefficient 
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Figure 3.45. 


Variations of the Reduced Pitching Moment 
Coefficient at Zero Thrust with ysina 
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Figure J.69. 


Left Rotor Pitching Moment Coefficient Versus 
Collective Pitch. Im « 90* Hover. 
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3.103. Left Rotor Normal Force Coefficient Vs. Rotor RPM. 
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Left Rotor Collective ^Degrees. « 15^. 
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Figure 3.131. Left Rotor Thrust Coefficient Versus Rotor RPM. 

IjT » 15^. Full Scale Airspeed 180 Knots. 
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APPENDIX A. APPLICATION OF LINEAR REGRESSION TECHNIQUES 


As mentioned in the Introduction < an attempt was made to 
apply multivariable linear regresi.ion techniques to the data. 
This was motivated by the need to process the large amount 
of data available in the shortest time A brief description 
of the procedure and the results is given heve. 


The basis of the approach was that each component of force or 
moment was functionally of the form 


Cp ■ CpQ (a, u, 0/ e 


where 


!fp, !fp 

»A, 3B. 


were 


‘1 ’"1 
independent variables. 


3Cp 

.75' ^1' ®1^ ^1 

, in turn, functions of 

i . e . , 



the 


(A-1) 



(o, M, ft, ^1' ^1^ 


(A-2) 


Thus each component was to be represented as the sum of the 
force or moment at zero cyclic, and perturbations in that 
force produced by the application of the cyclic controls. In 
principle it appeared possible to read all the wind tunnel 
data into a data file, specify likely parameters for correla- 
tion (for example, ucosa, uC^) and use the linear regression 
option of a statistical package CSTATPAK) to determine how 
well a linear combination of the selected parameters fitted 
the data. 


The first step in this process was obviously to use the linear 
regression feature of STATPAK to find the values of the cyclic 
derivatives dCp/3A]^, 3Cp/3B]^ since this would be faster than 
calculating the slopes ^rom the >;ind tunnel data plots. 
Accordingly those runs in which u/ a>,3 75, and RPM were fixed 
and A}^ or B^ varied, were processed by* the statistical pack- 
age. The equation for a particular force coefficient waf. 
therefore 


Cp ■ a b Aj^ 

where the slope, b, and the intercept, a, were given by the 
linear regression. The slopes were then plotted against 
ucosa. An inspection of these plots shewed considerable 
scatter in some areas. The wind tunnel data was reviewed and 
it was then observed that small but significant changes in the 
va3ues of cyclic and collective (that were nominally fixed) 
had occurred in a number of the runs. In view of these devi- 
ations it was decided to combine the Aj^ and B^^ sweep data and 
introduce thrust coefficient as a variable in the regression 
equations. The equation form chosen was 


A-1 



(A-?/ 


CF » CFq + CF^ + CFjA^ + CF^B^ + CF^A^C^ 


Using this form in the regression analysis resulted in extreme- 
ly good correlation being achieved. Figure A-1 presents 
comparison of sideforce coefficient computed using equation 
(A- 3) with the measured sideforce, as produced by variations 
in A and B cyclic at u = .1 and a » 80*. Agreement is to 
within 5% of the measured values. Figure A-2 shows similar 
agreement measured and computed normal forces at w o' .314 and 
a » 42*. 

Since good correlation was achieved using equation A- 3 it was 
anticipated that those portions of the final rotor math model 
equations which were to represent changes due to cyclic pitch 
would be written 

3Cf 

^ . CF. ^ CF, CT 

alf “ CF 3 ^ CF 5 CT 

The next step was thi n to determine how the coefficients CF 2 
through CF^ varied with u, a, u, 0.75 and hence obtain general 
equations tor 3Cjjp/3Aj^, aCj^p/sB^, etc., covering al3 flight 
conditions. It was at this point that difficulties were en- 
countered in find'lng suitable combinations of parameters that 
would achieve thri,. The coefficients could not be represented 
by a simple function of vcosa as is shown by Figure A-3 which 
presents the variation of the coefficients in the normal force 
derivative. Combinations of likely parameters such as usir. 1 , 
li sin a, yCm etc. were tried to no avail and it was the re fc re 
concluded that some theoretical guidance was required. Th. s 
led to the d<rivelopment of the analysis described in Appendix B. 


A-2 
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As described in Section 1.0 an analysis was developed to help 
provide some guidance in choosing basic parameters for use in 
the regression approach. Later the analysis was expanded in 
the expectation that the analytical estimates could be simply 
related to the measured data and hence provide a means for 
predicting the rotor behavior. While this approach met with 
only limited success, the analysis and the results are pre- 
sented here and may be used as a starting point for a more 
refined treatment. The objective was to obtain some closed- 
form expressions for the rotor forces and moments that could 
be rapidly evaluated by a computer. Closed form was desired 
since the resulting expressions could be inspected to see the 
significance of each term and identify simplifications. The 
analysis is developed based upon the following assumptions: 

(1) uniform induced velocity across the disc 

(2) linear twist 

(3) constant blade chord 

(4' first flap bending mode only 

(5) AO in-plane motion considered 

(6) no elastic torsion 

(7) linear section aerodynamics 

Because the analysis must apply to a highly-twisted prop-rotor 
at large axial advance ratio, the conventional small-angle 
approximation to the inflow angle at a blade section cannot 
be made. This results in a very lengthy development of the 
analysis and the' resulting equations are also lengthy. The 
treatment is an extension of that given by Dommasch in 
Reference 4 and makes use of some results presented in 
Reference 5. The final equations are presented in full with 
no simplifications made other than those stated. An assess- 
ment of the relative importance of each term was postponed. 

The analysis is developed in rotor wind axes i.e., at zero 
sideslip and at a resultant angle of attack related to the 
wind tunnel pitch and yaw angles by 

-1 

QL^ » Cos (Cosoi^pCosi;;^) . 


B-1 


'''' 



COMPANY 


NUMBER D210-11505-1 
REV LTR 


ANALYSIS 


Referring to Figure B.l let the deflected blade shape be 
represented by 

rc» i- 0) 

where S(x) is the first flap bending mode shape and 
/(f) is a first harmonic variation with azimuth i.e. 

2* rce t iLSi*) ia» -A, eoi>^ ~ (?■) 

Taking moments about the blade root, the centrifugal, 
inertial and gravity forces acting on the element of 
mass dm yield 

,i 1 4 Mi rCi &) 

where i»lk is the mass per unit length 

and O^&is the shaft angle measured from the horizontal. 

The aerodynamic contribution to the root moment is 
approximated by 

^ • t (L I 

The total externally-applied moment is therefore 

M» - Mi- M*r 

and the existing elastic moment is 

• • • '* 

L L Ifzc 

Now 

(<Mi» r2fUf‘ 

iHif frJlY'A-* 


where 


f, e a, 1 
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FIGURE B-1. Definition of Axes Systems and Quantities used in the Analysis, 
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I.’ ; I, * 


Mw ® [ * Mt/ ^Kff 


(i>\ 


where M^ is the blade weight moment. 


The aerodynamic contribution to the moment is now 
required. With reference to Figure B.l, assximing 
small angles for the induced inflow, the tangential 
and perpendicular components of velocity are 


Op • \r^U-C«ii' + ^/JIA 


(t) 


where 

and 

With 


SLU 
s V/Ml 

ij s Attiuf -a% fesf s' = ds/df- 


these components become 


Of * A + {Ca-*jS']jt.Coif t Sjj 
Oj * + s 




rORM 4«lt4 t2/«4) 
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The elemental thrust is 




duJL 


(it* iL a>s4> 

Ll * iP I// 4 cetg^dr 
dJ> * U Vf pqdr 
<ii * $-4 ■& a*^(9-4) 


L<0 

•4* d^p/tiU 


Jff 


Therefore we have 


dT ' Kt / Vf/Ui^ ~ 0^ ~ ^ 


K(^ Kf s laxc-iL^K^ 

4 * a/4 




This form for the thrust allows the high inflow conditions ot the 
tilting rotor to be treated in a straightforward manner. When 

thrust force the drag contribution to the thrust 
will be retained since in high speed cruise this can be a signifi- 
cant portion of the total. For the present purpose of calculating 
the aerodynamic moment, however, the drag will be neglected. 

The local value of blade angle^S^is 

P * Pe - fir ^ - /I/ to*4' -3/ 00 

: ff>. Ad 

S' * $» • 

^S * 4, f d, iMUtj/ 

4/ Otid S/ are the lateral and longitudinal values of 
cyclic, ^ the root blade angle, and is the net linear twist. 


MJUt. 
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Thus 


■& &0ets0' (n.) 

Cti$ -t US&' L9 aM.0* Of) 

for normal values of cyclic pitch. 

The aerodynamic moment about the blade root at aximuth f is, therefore, 

Ml : f 

+ x’A + f.*- p Ceftf' Cctrfs') i- Oi 

Wh6n performing 'the in'tegre'tion^integrels of the form 

j ^ S ^ ^ S n* ccseUi 

are encountered. The following notation is introduced to represent 
these twist integrals* 


PcnM 4 tiM 
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Nw 



Wm. , jrt 

» ( Ji* 

t" - 

r ^ ^ I 

T*^ * (j 

^ f * * 

• Coiffa 


® ^^69 J 




ii^UM / : J 

y * J * 2 


6 ^') 


^<0 


The values of 


and twist. 


I I 


are constants for a given mode shape 


with this notation equation (14) may be integrated. The resulting 
expression when simplified to include only steady and first harmonic 
terms is 

J!* • ^ T’jT -M 2 ^ -^/*A r/J -/I -C ^ 


FOMM 4 MM (a/« 4 ) 
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* -*<-7^ * tf tZ 

f A TJ tJS/ * \ -I^OtA, 7^, 

• 7^1 j C 

' - J (,7) 
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Equating this expression fci' the aerodynamic moment to the moments 
due to inertia, centrifugal force, weight and elastic bending it 
is possible to obtain the coefficients aQ, a^ and 9^2 
first harmonic flapping by solving the following set of equations. 


4e 4 C/idi. * f2| 


where 


' 


do 4 ^VL Ag 4 CtJ Az, 

C^f 4 # 4 ^jx, 4 4554*.* 

r 4 iA4|T^f 

i/'-r /:V / 
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Ex pressions for the Rotor Forces and Moments 

Because of the very lengthy manipulations involved in obtaining 
expression for the rotor forces and moments, only the first few 

<3evelopment will be given and then the final result 
will be presented. 

Thrust 

Retaining the drag contribution, the thrust on an 

(10) and (11), and 

coIffJcJentr-*”* products of the flapping 

IT* Kt[ 

+ 4*1 iV'^ 5' xl* C»yi £»*«/' 

+ jifi f s' Ui<i' t- Ji Sfs)!- 

4 Zrjy^ j- s' co^f 't zs s ^ 

exprlL*ion^fcrSe lIltuTT'' """ 


I ■ -1 

+ Ai r A, 4*^ 1 r ^4if, 


-J ( /J.4^ 


J 


40 

cs 


^OAM 4t2«4 (2 /«4) 
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- "C ■* i rlav) -A, >(<<=> C 

- 4. ^ y T" - 0,.^ n“ ] 

This equation is solved in the usual manner by using 
momentum theory to obtain a value of X consistent with the 
thrust. 

Rotor Power 

At blade radius r the contribution of the element dr to the rotor 


power is 


^ im 

* ICrM [ Of Ur Oj 

Performing the same kind of substitutions, simplifications 
and integrations as were required to obtain C„, the final 
equation for rotor power coefficient is t' « 

. i U3,-a.A,)TZ] ^ 

Cm) 
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Rotor Pitching Moment and Yawing Moment 

In Reference 5 a particularly simple set of equations are 
obtained for the pitching moment and yawing moment coeffi- 
cients in terms of the first harmonic flapping coefficients. 
Defining an analog of the Lock number as 


ft * (*5) 

the pitching moment coefficient is 

ar, 


and the yawing moment coefficient is 


CTlu 


r. 


where A| is the (nondimensional) first flap frequency and 
the directions for positive moments are defined in 

Figure B-1. 



Rotor Normal Force 

Referring to Figure B-1 the elementary contribution to the 
normal force at station r on a blade is 


With the approximation 

* di/M - 

the normal force coefficient may be obtained (after some 
considerable labor) as the sume of fi¥e terms 

iUJ r XCh9\ ^ iCyjlr 4 -f 

^ ^ ffX (TL 


FOf^M (2 
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Comparison of the Analytical Predictions with the Wind Tunnel 
Data 


At each wind tunnel test point, the test values of collective 
pitch, cyclic, shaft angle of attack, yaw angle, advance 
ratio and rpm were used in the analysis to calculate the rotor 
forces and moments. The first flap bending mode shape was 
computed using the model blade properties presented in Refer- 
ence 3. The frequencies of the first flap bending mode were 
obtained from the rotating blade frequency measurements given 
in that Reference. The average value of blade drag coeffi- 
cient, Cjj, was obtained by extrapolation of the hover power 
coefficient to zero thrust and using the approximation 

^d “ <36) 

The results are presented in Figures B-2 through B-6 which 
show the calculated and measured forces and moments. The 
calculated values are plotted using the same symbol as the 
test value but marked by a flag. A guide to the test runs 
is provided in Tables B-1 and B-2. 

In general the predictions are in fairly good agreement with 
thetest data considering the simplifying assumptions that are 
made in the analysis. In cruise flight the predicted slope 
of sideforce coefficient with collective pitch is opposite 
to the test data. This may indicate the need for a lead-lag 
degree of freedom in the analysis. Further work is required 
in this area. It should be noted that the execution time of 
the program incorporating the analysis is short, about the 
same as the curve-fit equations. 

It was hoped that the analytical results would be linearly 
related to the measured data i.e. 


C_ « a + b Cp 

THEORY TEST 


(37) 


where a and b would be functions of effective advance ratio, 
ucosa. If this had been the case, then estimated values for 
the rotor forces could be obtained from 


TEST 


1 

S 


(C. 


-a ) 


THEORY 


(38) 


In other words the analysis would be used to predict the 
forces and moments, and equation (33) used to correct them to 
provide estimates. While linear relationships were found 
for thrust and power, the remaining quantities were not 
linearly related and the approach was abandoned in favor or 
direct curve fitting as described in Section 3.0. 


B-18 


D210-11505-1 





u 

a 



EmSS 


_ «.75.. 

|[^^l||||||||| 

90. 


memm 

59 

__0 

a 

-n 

VARY 

mm 

■m 

7.9 

42 

g 

ii/ 



90 



8.2 

43 


*^.75 



80 

1 


VARY 

45 

A 




1 

VARY J 

8.9 

.46 _ 

.U. 




VARY. 

1 


I 

__ _50. _ 


L. . 



02 



k 

.. 5L . 

Q. 

LT5 



1 


i 


VARY 

53 

_o 

RPM 



i 


VARY 

8.9 

54 

A 

a 

.226 

VARY 



13.5 

55 

0 

.j - 



65 



13.5 

38 

i 

0.75 



. 1 . 

1 

VARY 

59 


RPM 



t 

VARY 

13.5 

61 

_o 

a 



VARY 



16.0 

L .62 

❖ 

- X^- -1 



_49.0 



16.0 

il- 

s 

. L75.. 



... I_ 

f 

VARY 


V 

RPM 

VARY 

* 

VARY 

16.0 

68 


a 

.315 

VARY 



22.6 

69 

0 





42.0 



22.6 

72 

_ o 

auem 



4 



VARY 

.7.4. 





i 

WBimm 

WBSKM 

15 



. a 

■■ 


mmm 



■iiKA 

76 


^ ... .. 

■■ 


_ 59.. 



18.2 

79 


0.7JL 



1 _ . 

f 

VARY 

ai 



mmm 



VARY 

18.2 

82 

A 

a 

.226 




12.10 

83 





76.0 



12.10 

86 


0.75 

i 


1 

f 

mESsSm 

88 


RPM 

VARY 

tf 

VARY 

12.10 



a 


VARY 

1065 

20.7 

! 90 


. >1' _ 

mm 


26 



m^mm 

93 


t^.75 

mm 

B 

. 1 

J 


VARY 1 

95 

- 


VARY 


VARY 

■EEKfl 

16 


JX . J 

.351 

VARY 



29.0 

97 __ 

A 








29.0 

' 100 

_* 

e.75 







r 

VARY 

, T02 



VARY 



VARY 

29.0 

1 111 

V 

a 

tmsm 

VARY 



32.0 

i 112 





24.0 



32.0 

i 115 


0.75 




__J 



VARY 

msam 



VARY 

l■U■i 


32.0 

.118 

iniii^Tjmiiiii 


.510. 


m*mim 

■Rffiil 

.119 


wmSm 

1 __ 





36.0 

moFfm 

€ 

mmm 

f 





l6.0 

Hi£n 


RPM 

VARY 




VARY 

37.5 

128 


a 

.58 

VARY 

830 

41.5 

132 

A 


.58 

0 





138 


a 

.707 



mm 

imm 

47.2 

141 

_ b 

. * 

.707 



■■ 

IHH 


.145 


. a . 

_ .836 







148 



.836 

■i 

jmi 

laimi 

mm 


153_ 

A 

a 

_J7 

imKiHi 

mEFFH 


154 

i _ 

. . <lL 

..9.7 




mcrm 

15 9_ 

1 

CL 

..4.48 




33.8 



- 

ik 

- ■Liig- 



V 

34.2 


TABLE B-1. List of Runs in which Cyclic was Fixed 
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RUN # 

SYMBOL 

RUN 

TYPE 

M 

a 

RPM 

®.75 

. 29 

__o _ 

Jil . 

_ .0 

. iO.O 

■IVIH 

■mm 

31 

□ 

B1 

_D _ . 

90.0 


a. a 

.32 

o 

B1 

__._0 . 

90.0 


11.8 

40 

A 

.Bl 

^1 

KIIPRH 

mm 

7.9 

41 _ 



A1 

.1. 

79.3 

mi:|cH 

8.0 

. 84 



A1 

_ .22 

76.0 

.1185 

12.0 

85 

. _a - 

.jar 

.23 

76.0 

nafi 

n .Q 

46 . 

.. o_ 

Bl ... 

HKHSi 

82.2 

■I'lim 

._8.6 

AQ 

2 


1 


mrm 

9.0 

56 


T~ 

. A1 

.21. 

mm 

mim 

■mmi 

57 


M 

Bl 

.23 

63.5 

1185 

13.4 

77 



A1 

.32 

59.5 

1184 

18.0 

78 


61 

.31 

60.0 

1185 

17.9 

64 


A1 

.23 

48.6 

• 1185 

16.0 

65 

— S— 

Bl 

.23 

48.6 

1135 

16.3 

70 


A1 

.31 

42.0 

mim 

22.5 

71 


Bl 

.31 

41.8 

mim 

i^mi 

91 

□ 

A1 

.25 

31.0 

1065 

20.7 

92 

■ 6 

Bl 

.25 

31.0 

1065 

20.8 

98 

0 

A1 

.35 

26.5 

mim 

21. 0 

99 

HWHH 

Bl 

.35 

26.5 

1065 

mm 

113 


A1 

.45 

24.3 

1065 

32.0 

114 

■PiflH 

Bl 

.45 

24.3 

1066 

31.9 

120 

HD9HI 

A1 

.51 

12.8 

943 

37.7 

121 

WKFiWM 

Bl 

.51 

12.7 

944 

■nm 

122 


Bl 

.51 

12.9 

945 

■am 

163 

S 

A1 

^45 

14.0 

830 


164 

■■i 

ii'HI 

Bl 

. 44 

13.7 

830 

■nraj 

133 



Jil 

.58 

15.0 

830 

■nm 

134 

■IS 

liHl 

Bl 

.5B 

mm 

830 

41.7 1 

142 

■H 

iH 

A1 

.71 

14.5 

830 

■nm 

143 



Bl 

1 .71 

■IHpH 

■mm 

■nm 

149 

HkOiH 

A1 

.84 

■mm 

^0 

■mm 

150 

■Bm 

Bl 

.84 

mm 

828 

51.9 1 



lAl 

.97 

14.4 

827 

■mm 

wpomm 


Bl 

..ai... 

1 14.7 

>?26 

mmi 


TABLE B-2. List of Runs in which Cyclic was Varied 
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CURE B -2. Thrust Coefficient versus Collective Comparison 
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FIGURE B-3 


.8 
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Predicted and Measured Thrust Power Relation- 
ship during Shaft Angle Sweeps. 
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APPENDIX C 


D210-11505-1 

LISTING OF THE ROTOR MATH MODEL SUBROUTINE 


1. 

\: 


I 


i 


j 


i 


i 

i 


Sl;fc?.CLTir.E; nCT.:.K(/Lt.TSXa;lTSJ*a:LL*ALM;K,H£7;.H»Rr 
IC'F»CSF*CP«) 
cixcNsrcN x(e«5*fe) 

C**** CALC EStiyATEC -VALUfci-OF FDRCrS A^D- MDr-EAt-T-f,- **.**.«t*.* ***•»- . 

LC 343 I=lf« 

Cn 333 J=X*? 

^ . — 

333 X ( I«u«K) :0»0 

X t l*-l - 

V (l*l*3)s-P74.53hl 
r<l. l*A)=ll99,6?fcf> 

X (U1 *.5.) =.-4-752^22 - - - -- 

X (1, lf6) =SS6.571fc 
ya*l*7)=-l£.5«1875 

X <2*^l*.l»-=--*0if.3- 

X <2. 1 *?) =1.‘>95S 
y <?. 1*3>=-21.6829 

X(2.-l*.<a=2.C,BAA.3 — . - 

XC3.1»1>=-1.3713 
X(3« 1*2>=16.9001 

- - XC3*l,.3»i^.2La-i'i&£ 

y(3*l*A)s'17,t2Sl 

X(3*l*‘>>=-f03.7611 

X<3*.1*.&).=-177.^2S^ - 

XM.i»l>=.7213 
X(<». 1.2) =10.«Alfc 

X(*»l,4)=6flh.3?Al 

X(4,l,!;.)a-823.16aH 

X <4*1 

X(5*l*l)5-.926l 
y<5«l. 2)2-15.6913 

- X tS.»L»3).s:U6.>06S.7 

X('..l. 4)5-150.7721 
X(5. 1,5) =t0.ft424 

. X(6,.l.l-)-2i»£.0a2 — 

xr&.l«2)2-.7260 
X (6,1 .3) 566.1302 

. Xt6,.1^4>JS.il^9*5Ai& 

x(f.l,5,)=203.539H 
XCt.l, 6)5-70.6067 


- oucncoic 

OOOOCC2C 
00000030 
■ 00000040 

0CQ0005C 
OUOCOOhO 
00-000070 
OCOODOeO 
OOC00090 
&0000100 
OOCCOllO 
OOC0012C 

- 00000130 
00000140 
0000C150 

OO0C016O- 

OO0CO17O 

ooocoieo 

OOOC0190 
00000200 
OOC00210 
-00000220 
00000230 
00000240 
00000250 
00000260 
30000270 
--00000280 
00000290 
00000300 
- 00000310 
00000320 
30000330 
.00000340 
00C00350 
0000036C 
00OC0370 
00000380 
00C00.XQ0 

-- -00-000400 
00000410 
00000420 


PAoe » 


C-1 



i m 


D210-11505rl 


X<7,I*1)=21.6L21 . - 
i' (7, 1 ,?) =gC*9fc'^6 
:M7»1.?)=-I0.2e7b 
.<{ 7 , ’.^•}s.-»77.49.2C_-. 

=fc,7245 

. 

X<8.1«2)=-I41,76fc7 

X<ttl*4)=/>67.<5147 

U.L*-l*^X=.rLlZ7..t7J>J. 

X ( =-.1416 
X(1.2.2)=-f^7,2749 

- X^(-I»-2-*Z.l-=2<U-.9fc.59 

X<1.2*4) =-46C.43tJl 
XC 1«2«5) =399.693<i 

X( — 

X (2,2 ♦!) =.G267 

— . X (2,2,21 =-!'2*tS25- 

V (2.2,3) =.f.r.50 
X(2, 2, 4) =-14.6560 

— xta.-2.,.5iJu;L-x.as& 

X(3,2,l)=-1.4225 
X (3,2,2) =5.1466 

)-= *2-7 .-420-1- 

X(4.2,l)=-1.2323 

X(4,2,2)=1.7646 

X-t-4,.2..04-=.-.46.»56.(iA 

X (4,2,4) =63.0376 
X (4,2,5) =-28.4635 

Xt -)-= = 3-.X12-4- 

X(5,?,2)=-5.1175 
X(5,2, 3) =-11.7354 

)US,-2..v4*x40»65« 

X(5,2, 5) =-26.7432 
X(6,2,l)=-.9117 

<-6-r^^P )-=...t7..44-9^ 

X(6, 2, 3) =106.6433 
X(6, 2,4) =-122.5628 

) -=4 0 > 0 6>^ 

X(7,2,l)=22.9366 

X(7,2,2)=43.6971 

5t.<Z,^»3) =-261L*.a507— — 

X(7,2,4) =-257.9979 
X (8.241)=3.3847 

X-(6^,4J.^S.l-=4-OA-.OOa3 

X(8,2«3) =-854.6753 
y(8,2,4) =2363.6416 

X-( a,-2-^Sl = -,24-50*.8-9 1 4 

X(8,2,6)=907,4531 
X(l, 3,1) =-.0711 

X4-W3-.^-1-=.1-4S - 2 - 74 1 

X(l,3,3)=-1732.9126 


CCCC0430 

00 000440 
0000D459 

— OCO-OC460 

0000047C 
00000480 
oacoc49c 

00000500 
0OCOO5 10 

CCC(V&52 0 

00000530 

00000540 

COCOC550- 

DC00056C 
C0000570 
— —— — 0 do 0-C 5*8 0* 
00000590 

0OO0C60O- 

00000610 

QQ000620 

00-0410630 

00000640- 

00000650 

00-00006-0 

00000670 

00000680 

; — 00-000690 

90000700 

00000710 

000007-2 0- 

00000730 

00000740 

-00-000-7-50 

00000760 

00C00770 

00 00-078 0 

00000790 

00000800 

00000810- 

QQC00820 

00000830 

00000840 

00000850 

00000860 

00080870 

OOOC088C 

00000890 

000089-0 0- 

00000910 

00000920 

00-000-930- 

00000940 






i 

«v 




t 


C-2 


V. 
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XCl*3.r) =-?1977.79:l 
Ml* ‘«i.)s2f;^02.0fe95 
X <1 *3*7> =-l‘:9i!9*1633 

i‘(l*»3-*-^)-*75.2^4 70 

X<2*3*1>=.0110 

X<2*3*2>=l.fe510 

X(2*3*^>sP.50b7 — 

V(2*3*n)s3.e*l397 

X(2,3*r.)s-220.2950 

X(?*^3*(v)-s2bT.-7575 

X<2*3*7)s-115.7fe26 

X(3,3*l)=.2fc53 

X(3*.3*2'JS.-7^5021- 

X(3*3*3)s34.26b7 

X(9*3,l)=-.0ll7 

. 

X<**3*3)=78.412fc 

X<<‘.3*4)s-lb7.1991 

X»4*3*e-X-74-S577-- - 

X«^*3*1)=.3494 

X(5*3*2)r4.5443 

- -X<-6*-3*.3)s.-7»2 0Sl« 

X('^*^*l) = *0b79 

y(f.*3*.2)=b*2§.19 

X(6*3*?)s-56*6109 

X(6*3*4)sb7.6128 

. X t f *-3-*-5 

Xt7,3*l)s2.l864 

Xt7,3*2)=40.1364 

- X(b*.3^1)-=-l,^9i4. 

X(8*3*2)=39*4773 
X (8*3*3)5-308.0650 

X (H.3*Al-846»..^53 

M8*3*5)=-978.*8486 
M8*3*(‘) =44b*9769 

X.(3*-4*.U-5.0208 

X (3*4*?> 52*6330 
X(3*4, 3)5-28*2534 

X(3*.4^7s49.AJ62 

X(3***5)5-11,0460 

X<4,4,1)s.1416 

X (4.4 ♦^Xi-5♦042.6 

X(4, 4*3)533. 9888 
X(4,4,i»)s-1C*3021 

X(5.4.tXJt,4Si5 

y<:*4,2) 5.2,9718 
X(5*4, 3)529.7889 

X(5*4*^Xi.^..^l4 

X(5*4*fi)=l0.4779 
y(6*4*l)5,l063 
(6.4*J)£5..JLiA2 


noo'oC-5? 
OQOOO^fcC 
0CC0C97C 
000009HO 
- - 00000990 
000010C3 
00001C13 
000C1020 
00001030 
0OCO1C40 
OC001050 
00001060 
OCOOIC7C 
OGG01080 
00001090 
OOOOllGO 

- (WOOlllO 
00001120 
CC001'13S 

• 00001140 
ODOOllbC 
00001160 
00G01170 
00C01180 
00001190 
0300120G 
00001210 
00001220 
00001230 
0C001240 
OOC012bO 
00001260 
00001270 

00001280 

0CC0129C 

OC001300 

- • -OOC01310 

00001320 

00001330 

-■-0000134C 

00001350 

90001360 

- 00-0015 70 
30001380 
00001390 

- 0CC014j0 
OCCC1410 
00001<^2C 

- 00001430 
00C01940 
00C0145C 
■ -00001460 
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y (f ♦«»3) =-12.0f^5? 

-X < 7* ^ ♦ 1 ) =»03fc0- - — 

y(7,*.2)=<.e.3hlt 

X(7*4.3)=-93.3694 

X(7*A.4)=3S3..3Il£i - 

X(7.<f*f )=-6?h.lf!5 3 
X (7,A*6)=293.,?343 

Xlh»4,l)=-.C303 

X(fc.4,Z)=6,f2^1 
y<f *«*3>s-23.y^3 74 

Xt£*-^*.^J.= £C*£fc3fc. 

y S-3£.17C2 
> ( I* 5*2) = .94 

X.<2..£..ll=-:i.C.C.aS7 — 

X(2.5.2)=-.007»4 

X(2.5.3)=C.89403 

Xti*3iU).=-aCJ}-lt 

y (3.5.2) = .C234? 
y(3.5«3)=-l. 12935 

.Xt3«^b.^4.Xsl..71.22i 

X(3. 5.5)=-. 74232 
X ( 4,5.1 >=.OCC 55 

-jLt-4-,ii.,^l^--.tU^lS 

y (4.5,3)=-5.2fi9fi3 
y{4,«^,4)=4£.19186 

,.X,(4,5-,.5)=J’-12£..Z31LZ4 

>(4,5.6)=141.12570 

X(4,5,7)=-55.68757 

. X ( 5 * s.»x.)Ls.=.,.Qa.a.i.a 

X(5.5,2) =-,09926 

Xt5 ♦ 5 , 3 ).= . i2A X4—_ " 1 — 

X(5,5,4)=-.C5645 
X(t,5,l) =.0005283 

_Xtb..5.2> =..1-723*45-9 

X(6, 5, 3)=-. 7342283 
X(6,5,4)=l. 6377934 

.X<6.5-.5)---l .-7-335-0-4^ 

X(6, 5,6)=. 6787662 

* NOW FORM DERIVS URT A 1 *8 1 ♦ ALF A tPSl ,£TC ****** 

-Z-AS-S ( AHUF^C-O-S OL ’fiHR-) 

WsAMUF- 

CTAlsUK(X,Z.l»l> 

CFA1=C.KAX*242»-3> 

C\FA 1=GK (X*2»3»l ) 

CSFA1=GK(X.Z,4,1) 

CP M A-liilX 4 X 

CYFAl=r.K (X*2*6,l ) 

CPXA1=0K (X,2,7,1) 

F b AX .2 .-f ,^1-X 

CTE1=0K(X,2,1,2) 

CPElsQK (X,2»2,2> 


00001470 

00001480 

00 001*9 0 

00 00 15 00 
00001510 

— occcir>20 

00001530 

0000154C 

00CC1550 

OCOOlC’bO 

00001570 

- .-00(U!15«C 

00001590 

00001600 

oo-&oi6-ia, 

00001620' 

0000163C 

4)20016-4 0 

00001650 

00001660 

— 00001670 

00001680 

0000169C 

04)201700. 

00001710 

00001720 

-02001730 

00001740 

000C1750 

2 22 23-7 6 2- 

00001770 

— :: -00 001 780 ■ 

0CCC179C 

000018CC 

20001810 

00001820 

00001630 

O02018A0 

C0001850 

0000186C 

00001-872 

OOOOlHhC 

000C1850 

20021900- 

000C1910 

00001920 

ao 2219 30 

00001940 

00001950 

02001960 

00001970 

00001980 
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C\F"l=r,K <x»Z»3»2) 

CSFr-i-c. (»,Z*4,2) 

c?' o'.z^etS) 

CYV- Is'-K <X*/ ♦C,2) ... 

C?-I-1 = CK (X,Z*7f2) 

Ftr^'rlsOK (X*Zffct2) 

CTALFsr.it — 

CPALF=r.K <x*2*2*3) 

C^FALF^r K (X,2»?»3) 

C'«PALF=AhS(C*\FALF> — - 

CSFALF=GK(XtZ*4f 3) 

Cr’'ALF=f.K(X*2*Sf 3) 

CVVALFs<iK(X»Z^<^*.3) ~ — 

Cr-''ALF=f K(XtZ«7f 3) 

Fr'AL^sGK(X,2»B»3) 

CTF3lsCtc<X»2f-lt^) - - -- - - 

CFFSTsGK (X.Z«2«4) 

C^FPSI=^.K (X,Z*3*4) 

CSFPSI=CKtX»Z*-4.,^>— . - - — — 

CP‘'PSI=0K<y,2»S,4) 

CY''f’SI=GK(XtZ*6*4) 

CH ' f SI=ttK<X*Z»7»t^>. 

FPvrslsRKfy^Z^e*^) 

CTcpvrCK (X,Z«1»5) 

CPF?K=GK <X*.Z*-2*&> - 

Cr.rPP»*5GKCXtZ*5*5) 

CSFRpysGK<Xt2*4,5) 

CF.VPPK=GK<X»Z*5*^5) -- - - - 

CYyRPysGKCXtZte^S) 

CE'^RPM=GK<X*Z*7*5) 

FBvprMs&K-AX«-2«^9.&) 

%CW CALCULATE REFERENCE VALUES CF ALPHA ***•*•• 

ALF=:EF = F9,7l76*Z*(-322.e 042-»Z • < 77C .9 1C«>Z * ( -1 A * 6 . 3424^Z« 

1 (lAA9,9&ij^4.Z*(-54.5,ll54) >)^)> 

SKLFACsl .0 
hpvrRFP.SKLFAC 

£YE^•ACsO-*.^- - 

ccys-tFsl lfcE.-<llfe5.-fi30.)*( l.-EYENAC/*'.) 

IF<tYENAC*CT*45*-> RPKREFiUt.5^ 

Cf LF-PPsFpP-RFfIBEF 
CELAlsAlTST-5.0 

UELblrSei.WT-'tS*^ 

pelalfsalphr/ctr-alfref 

SELfSIs-EETAR/OTP 

APSIsAPSCCELALF-V 

AfSSsAFSIOELPSl) 

rELCT = (CTAl*CELAlt.CTRl*DELBl^CTALF*CELALFt-CTPSI*AeS2A.CTRPM* 

-1 CELRPKXA4>aft&0.». 

CELCP=<CFA1*0ELA1^CPB1 *nELBl*CPALF*D£LALF*CPPSI*AE«S?*CPRPH 
1*CEL»FP)/X0000. 

CELCNFs(CNFA.1*&ELAIaCNFE l*l>rLf'lAC.'.FALF*CELALFt.Cf FP^S-I *0ELPS1 
1*CNFRPF*PELFPF)/10000* 


00001990 
OCOCZDoO 
COOOZ'JIG 
00CU2C2C 
00CC2030 
90002C4C 
• OOOOZCjC 
ooooauoo 
00CC2C7C 
00C02060 
OOOC209C 
00G02100 
- - 00002110 
00002120 
00002130 
- 00002140 
000021BO 
0000216C 
00002170 
00002180 
00002190 
00002200 
00002213 
00032220 
0000223C 
30002240 
00002250 
00002260 
000Q2270 
30032280 
- - 00002290 

00002300 
000C2313 
00002320 
C0C02330 
OCOC2340 
00002350 
03C02363 
00002370 
000023<'0 
00002390 
- - -0000240C 
O0C02A1C 
0000242C 
00002430 
00002440 
000024b0 
o-ao&24(i,c 
0C00247V 
00002A(r0 

00002A90 

0000250(1 


C-5 
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[;.tLCSF = (C.':‘^Al*0C:LAl4C£.Ft l*D!TLf-l + CrrALF*:' L.‘Lf->»CSFPSl »DuLPST 

l+CSrRp:'*rFLFF?n/lCftOO»- • - - — 

i;FLrpy=( cM'Ai »:^E:LAi>crr'Hi*DE:Lr-i*cr f--AL? ♦of.lalp^cf f'Pr.i^LiLPSi 
l^CF? RP.‘-«DtLFPf- )/lCD00. 


00C02S1 0 
OaP025CO 
00C02^30 
00002540 


CELCY>"=( CYJ'Al*r.£LAl>CY:<Bl*DELFl*CYPALF*C.ELALF-»-Cffc!F-Sl *DCLPSI - ■ - -0CC02fib0 


l4CY!'So^*r.ELFFf--)/lCC00. 

DELCor-<'KFAl*CELAl*CPMBl*0ELBl*CF"ALF*ArSl-»C6MF?iI*ABS2 

.. . -1 ♦CL-.*'RFK*LF.LKPt’ 

DELF^''=F?f-Al*CELAl>FEf'*?l*DEL81^FPPALF*ABSl^FhVPSI*AFS2 
1 ♦F“’‘«’Pt-*GELPPF 

IF(7,LT.C.G15)~G0 -TC. 30 

IF(/,LT.C*C55) GC TO 31 
IF(2.LT.C.100) GO TO 32 



IF(2.LT.C.lfe2) GO TO 34 
IF<2.LT.0.210 ) GO TO 35 

- IFA2*L.T.O*23C)-.&0.-tC-3G- 

IF(2.LT.C.234) 60 TO 37 
IF(2.LT. 0.410) 60 TO 3H 

IF(2.L-UG*4SJ)._Wk-Xa~I3 

60 TO 40 

3C GRAri=11.533 

. - - CEPT1 = .4'.£2 

GP Ar2=106 .FG7 
CErT2=-l.f 

- GF.A.a3.=.f.2£7 

CErT3=.054 
GRA04=-G.f 67 

-4tr-T4=»0-i5 

GC TO 41 

31 GFAC1=3.175 

CEP-U.5 L.IU2 

GRA02=92.75 

CEPT2=-1.291 


00002560 

00002570 

C-0C025PC 

00C0259C 
C0002e00 
• — 00002610 
0C002620 
00002630 

00-G0264 0 

00002650 

0C002660 

-000026-70 

00002660 

Q0C02690 

OC 0027 00 

0000271C 

C0C0272G 

00002730 

0000274C 

00002750 

-00002760 

CCQ02770 

00002760 

— : 00002790 

00002600 

00002610 

00002620- 

00002830 

00002640 


- -GRAD3R..Z2L5-. 

CEPT3=.C87 
6R AC 4=5.5 


— . .. -CEi.T-4.;Ra...I.9I 
GC TO 41 

32 GPACl=1.467 

CEE LUO. UU 


GPAD2=77.556 

CEPT22-.456 

— 66mG 3 s. « 5. 7 6 . - . 

CEFT3S.C67 

GiaCL4.s.l.-l-l-l_ 

CEFT4S.C49 
GC TO 41 

13 G-RAClS-«^I6X_ 

CEPTU1.333 


GPA:2s65.106 
C6P.T2=.-l*2JU 


00002850- 

00002860 

00002870 

00-002880- 

00002890 

00002900 

00002910 

0000292C 

00002930 

0-0.0.02-8.40- 

00002950 

0-0-802-960 

00002970 

00002980 

08002-990 

00003000 
00003010 
8800 1 0 2 0 


C-6 


. -rf- 



GRA:.7 = . = r7 
CC‘ TTs.rs'; 

GS4C4 = 1»7C£ 
CfPT*r-.ClC 
GC TO t,] 

GRAClsiC*4 

CCFTlz-,312 
Gf<4L'2 = 193.I33 
CCFT2S-J 7»i2 -- 

3 = -. 733 
C£FT3=.?7h 

coAr^s-10. 

CEPT4=i.71 
GO TO 41 

_ 

CEPT1=1.768 
GSAD2=4T .633 

CEFT2=6.7T5 

GCAC321.47G 

CEPT3S-.081 

CR-AC-4J4?^1-2-5 

C£='T4s-l ,22£ 

GC TO 41 

GRAGlat^lM- 

CEFT1SU009 

GRAr2=-40. 

..CFB-T2=*A*^ 

GPAL3rl, 

CEFT3S.C20 

-- 

CEPT4S-.780 
60 TO *1 

- 6RAfi4 s-5S-^ 

CEFTi=-lc.H02 

CRAL2S1125. 

CEf^T2=.*2 4^.45 

C4AD3=-1C, 

C£PT3=2.05 

aR-AR-A-i-grOi^M 

CErT4s24.463 
GC TO 41 

6RAC12-W 

CEPTlsl.249 

6RAC22G1.8182 

■ CEPr2Jt5...83.5 

GAA03S1.364 

CEKT3S-.IC9 

. GRAO-A-r^ZfvZ 

CEPT4S.006 
GC TO 41 

GRAris3*^C6W>. - - 
CEFTla.ACSA 


59 


D210-11505-1 



f 


41 


21 


2 ? 


2i 

?4 




C210-11505-1 

7.7C1 . — 

C£* ■'■:r7.3l'2^- 

G^A24s,793 
C•^•'To = -.C^^ 

CC TC 41 - - 

2PAr l=l.h79^ 

CCPTlsl.0901 

.Itt- 

C* M?sl C ,‘:7‘^9 
GPiPZs.ftfe 

CETT ’ = »&«£._ - — 

pc AT4S.743 
C£'PT4 = -.005 

. CE AncTs-c GRAC i^r^cErmvi CCO.-. 

CEI- TCT = GP An2*Z^C£FT2 
Gc Ai:CP = GR AC2«2-»CEFT7 
c £ r T c r = t G r. ad. 4.* z-^c £ f t / i g o 
IF(W.GT.C.51) GO TO 21 

C£FC\F=-1.4/(95>V(*n, llf4 + W*(4fi,14f’.r + W*C46.1E3?'»y*(-41lP.74 72 
1 ♦'** (2711.4 •Su.92*W* (-1478G,65‘V9*W* (G9504,RCI-5*W.*4 -59 c‘18-»^JG83 - 
1 ) 1 ) ) ) ) ) ) 

C£FCr,F=CCFCf F/IOCC* 

-fcP CC\Fi-.,0O.‘-^7-tW-* < 2^12 2C-»V*^t-^l G.9 482*U« ( (4-2G^ 

tC14l + W*(-1 043.6fc43^V.*(C2 0.fc263*W*(lG2«.?7744£*(-g6G, 

1.7f 33) ) ) 1 ) ) ) ) 

..XC--ia -22— 

CEFC\F=< .75-35.79* < U-*S1 ) }/l 000 . 

GRCCr,'F = ,0 3Q46£7-.264 94*( W-.ri) 

- iFtk.G7^c-.x;:^5i— C7i«:..'.F-a. 

IF(W.GT.0.50) GO TO 23 

CFFC5F=-.31fe4*y*(31.83C6*W*(-457.3fc74*U*(4745.5645*W*(-3 

-166^a.2.UL7.'».WA_t4,la.a9.1^22h6*U*( -.44 V35.S..4a4fr-«Xj»-(.54i(U^4 

1 .9029^U* (-3C3701 .0«C9) )))))) ) 

GO TO 24 

. CEFC.SF=--*343.2.7.tli.*9*.<U-.5J 

CEPCGFsCFFCf-F/lOOC. 
lF(w,GT. 0.226) GO TO 25 
-G.RC CSE.s-!^iIlZ*ktjiLt.-X*6J:3.4.*W.4 (...0a&.lAfcU.(.7 &.< 

1) 

GC TO 27 

- 1 F tU .tl *1U554.-X10-XC~ 2 6 


4.<tia*.W.*X*.ULS^034 


26 

27 


GPCCSFsO,*w* C-12.830 6^W* ( 213 . 46454U* 1-11 72. 834f*g* ( 16 16. 760 2 
1*U*(4765 .9642*W*(-207fi2.560P^U*(26637.244r*W*4-11791 .4688 

.!_).) } ))-UJ 

GO TC 27 

6F(;rSF = .2119 56-i.3 5* (U-.55) 

( W..C.l..lJ.^7a7.1._SKDCSF5.0.. 


CEPCP”=-2.3F61*W*(34.23164W*(-266.fi538 4'. 
1 7636 *W* (-12450.4146 ♦W *(11267, 6934) ))) )) 

- -I K t y . G T *C ,..4.)_XrF 


*(2Sa.2259*W*(3745 


0 C C C 3 r 5 c 
0 f 0 C 3 5 6 r 
00C03570 
-- -OCCOZ-EftC 
0CC03?a0 
30CC36C0 
- -00003610 
0CCC362C 
0CC03630 
-■•&CC0.3t-40 
00CC365C 
OC00366 0 
— 0CG0367C 
CO CO 36 6 0 
00C03650 
0(»003700 
00003710 
OCCC3720 
00003730 
00003740 
00003750 

— 00003760 
CCQ0377C 
00003760 

— 0OC03790 
00003600 
30C0381C 
— aOC03F20 
OQ00363C 
00003840 

— --000038 50 

0C003H60 
00C0367C 
— 00003660 
00003690 
000039CC 
00003410 
00003920 
00003933 

— 10003940 
00003950 
00003960 
OC00397C 
000039KC 
00003990 

.80004000 
0000401 C 
0C004C20 
■ - 0CC04030 
00004C40 
00CC4C5C 
00004060- 


-i 


C-8 


( 
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ctr CP* iCf:i-Cf r-/icoo. 

IF(* .GT.n.ilt) G<»tJCPr-sO. 

IPC^.GT.C.75) 60 TO 2f> 

CKl:CY''=r •.CC-fc-S + i.^A tK 7‘>»r G — 

193^v •C-13C3r.fc<'fc2»W*(5919C.r 7C7*W*(-lJ'7927.£r-:t *U'*(10br<*5 , 

1 7'f.i )))>))>>- -- - 

Gf T'J 2't 

rp GFi’CYPslf .99P1*W*(-131 .5322^W*C3fcF,993 0^U» (-?t;4,3 72P^W* 

— 

?9 CCrCY-^sl •lfc2*U*(19.fc6«2^W(-lfe«.lP344.i,*(a??.12f c^W*(-3 

1 3. lt73>',.*<7979.29«l^W* <-lC0fc'’.b?9 7-**,« < .P>. U^u* 

1 <-177*.-bl29>K))^)^>')-) — - - 

CErCY«'sCEPCYI-/1000. 

CTtASEsGfiADCT#(COLL-CEPTCT) 

CFFyASE-sOK^&CP*G-T&A.St-4€gF-T^ 

CNFr AS=CRCC? :F*CTBASE*C£FCNF 
C£Ft AS=GFt CSF*CT£AS£^C£PCJ;F 

CF *'PA4;.s6.R.tCF-MAC-T€uA46-kCEPCPA:- 

CYvp AS^6FC'CY*‘•CT5ASe♦CEPCV« 

CT=CThASE♦DELCT 

- CPsCPRAivE ♦DCtrCP-— — •- 

CSF = CSFF AS*jtLCSF 
C^FrC^F^^ AS*CELCf»F 

Crf'-=CP»*-F-AS*4^tCP>» 

CYP = CYFPA£-*CELCYF» 

RETURN 

FU^CTIf^ OKCGfDUf**It J) 

Olf^ENSION C<6t5»«) 

GKiT. < I 

1 I,J*6)*rU>**Cf.( I,J,7)>CUH*0( I* J.fi) )) ) ) )> 

RETURN' 


C-9 


* 


3CCC9CT0 
0O0C*C80 
C0CC4C93 
3C OC^l OC 
C0C04 1 1 C 
OOC0412C 
CCCC^K D 
00GC414Q 
OCCOAlbC 
CC0C41fcC 

- OOOOVlTC 
CGCOAIFI 
0C0C<*1''’0 

- - 0000*-200 

D0C0A21C 
0000*223 
- ' ■Wr0&4-236 
C0CC*240 
0000*250 

0&004.260 

5000*270 
- 00 004280 

0000*290 
0000*300 
CC00431C 
-- 00 0 0*32 0 

0000*330 
0000*340 

S&C0*350 

OCOO*3f-0 
00004370 
0000*380 
000043 -n 
5000*400 
0000**10 



